Extracting Multidimensional Phase Space Topology from Periodic Orbits 
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We establish a hierarchical ordering of periodic orbits in a strongly coupled multidimensional 
Hamiltonian system. Phase space structures can be reconstructed quantitatively from the knowledge 
of periodic orbits alone. We illustrate our findings for the hydrogen atom in crossed electric and 
magnetic fields. 
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Periodic motion is widely recognized as the most 
prominent feature in a wide range of dynamical systems, 
from astronomy 0,0, 0, molecular vibrations 0, 01, 
chemical reactions [(j, particle accelerators 0, atomic 
and molecular physics |a |3, ll(J to fluid dynamics, e.g., 
statistics of turbulent flow Because periodic or- 

bits are the fundamental building blocks of the dynam- 
ics, they offer a road map to the intricate geometrical 
and dynamical structure in a multidimensional phase 
space. Unfortunately, this map is hard to decipher: High- 
dimensional systems usually possess enormous numbers 
of periodic orbits whose geometric appearance in config- 
uration space gives no useful hint at their systematics 
0, Il0l | . Indeed, it is not even clear a priori that a sys- 
tematic organization of periodic orbits should exist at 
all. 

The hydrogen atom in crossed electric and magnetic 
fields is typical in this respect. It serves as a paradigm 
of strongly coupled multidimensional systems because it 
can be investigated both experimentally and theoretically 
[UEil El However, after two decades of intense scrutiny, 
the overall phase space structure still defies a complete 




ol -" ' ■ ■ ■ ■ ■ ■ 

10 20 30 40 50 60 70 

period T 



FIG. 1: (color online). Periodic orbits at E = — 1.4a.u., 
F — 0.5 a.u. FPOs and their repetitions are shown with black 
circles, T% POs with red plus symbols, T% with green trian- 
gles, T£ with blue crosses, and T3 1 with magenta diamonds. 
The inset presents the 3-torus POs in an enlarged T-S range. 



understanding. In particular, no pertinent symbolic dy- 
namics are known. Although extensive lists of per iodic 
orbits (POs) are available in the literature |9L lici. IXia . fl4 | ■ 
no comprehensive ordering scheme for POs in three de- 
grees of freedom has been proposed so far, nor has it 
been shown how POs can provide insight into higher- 
dimensional phase space structures. That such a scheme 
should exist becomes clear from Fig. Q] We find that the 
periods and actions of most periodic orbits (POs) in the 
crossed-ficlds hydrogen atom fall into clearly discernible 
series. It is the central aim of this paper to elucidate the 
dynamical origin of this remarkable structure. We will 
show that a hierarchical ordering of periodic orbits arises 
from the hierarchy of invariant tori in an integrable lim- 
iting case, and we will develop the tools needed to char- 
acterize this organization in a non-integrable system be- 
yond two degrees of freedom. In the course of this work, 
the nomenclature used in Fig. [I] will be made clear. 

The Hamiltonian for the electronic motion in crossed- 
fields hydrogen reads, in atomic units, 

H = \v 2 - \ + 1 (p y x - Px y) + ^- (x 2 + y 2 ) - Fx. (1) 

Here r = (x, y, z) are the usual Cartesian coordinates, 
p = ipx,Py,Pz) the conjugate momenta, and r = 
y^x 2 + y 2 + z 2 . The magnetic field B points in the z di- 
rection and the electric field F in the x direction. Exploit- 
ing the scaling property of the classical Hamiltonian QJ, 
we can set B — 1 a.u. without loss of generality. We take 
F = 0.5 a.u. and consider two energies, E = —1.5 a.u. 
and E = —1.4 a.u., slightly below and slightly above the 
classical ionization threshold E\ = —2\/F. To identify 
POs, we employ a simple shooting algorithm, which like 
any other numerical PO search is not guaranteed to find 
all POs. However, our method of torus construction does 
not require the knowledge of the complete set of POs. 

To illustrate the mechanism that leads to the hierarchi- 
cal ordering of POs, we will first study a two-dimensional 
subsystem of the crossed-ficlds hydrogen atom, namely 
the x-y-plane perpendicular to the magnetic field. Its dy- 
namics can conveniently be discussed with the help of a 
Poincare surface of section such as Fig.El The plot shows 
two elliptic (stable) fundamental periodic orbits (FPO), 
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FIG. 2: (color online). Periodic orbits in the Poincare surface 
of section y = for the planar subsystem, E = —1.5 a.u., F = 
0.5 a.u. Non-fundamental POs are labeled by their winding 
ratios. 

called S + and S~ , surrounded by POs of larger periods. 
The latter arise in stable/unstable pairs according to the 
Poincare-Birkhoff theorem [l5|. (For graphical reasons, 
only one orbit of each pair is displayed in Fig. [5]) They 
can be thought of as being generated from the breakup of 
a resonant invariant torus in an integrable limiting case 
of the dynamics. In the situation of Fig. we say that 
the FPOs serve as the organizing centers for the longer 
orbits that surround them. 

Similar structures can be expected to exist in higher 
dimensions and have indeed been observed in integrable 
systems 0. Given that in a Hamiltonian system of L 
degrees of freedom there can be invariant tori of dimen- 
sion up to L, one might expect to find an entire hier- 
archy: FPOs serve as organizing centers for families of 
POs that arise from the breakup of 2-dimensional tori 
(2-torus POs), which in turn organize families that are 
generated from 3-dimcnsional tori (3-torus POs), and so 
on up to the maximum dimension. However, Poincare 
surface of section plots are hard to visualize beyond two 
degrees of freedom. Therefore, they do not provide a 
practical tool to identify these families of POs and to di- 
agnose their mutual relationships. We will here introduce 
a quantitative method to accomplish these two tasks in 
a multidimensional system and we will demonstrate in a 
3-dimensional example that the hierarchy outlined above 
does indeed exist. The hierarchy of POs that we derive 
mirrors the hierarchy and topology of their originating 
invariant tori in the integrable limit, which is also the hi- 
erarchy of surviving invariant Kolmogorov-Arnold-Moser 
(KAM) tori of the full non-integrable dynamics. 

The key conceptual tools that we use to elucidate this 
hierarchy are action variables and winding numbers. In 
the integrable limit, any motion is confined to an invari- 
ant torus that is characterized by a set of conserved 
action variables I. Their conjugate angles 6 determine 
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FIG. 3: (color online). The hierarchy of n-torus POs in the 
crossed-helds hydrogen atom at E = —1.5 a.u., F — 0.5 a.u. 
Insets show the FPOs in configuration space. 



the position on the torus. If an n-torus carries POs, in- 
teger winding number Wi , . . . , w n can be assigned to it 
so that the angle 9i runs from to 27r • w; before the 
orbit closes. Because winding numbers are topological 
properties, they remain meaningful in the non-integrable 
case, where the resonant tori have broken up into isolated 
POs. Using the winding numbers w, the total action can 
be written as S = w T(w), where I(w) are the action co- 
ordinates of the originating torus. In fact, the actions de- 
pend only on the ratios of winding numbers, e.g., w\jwi 
and W3/1V2 for 3-torus POs or w\/w2 for 2-torus POs. 

While the assignment of winding numbers to a given 
PO is simple with the help of a Poincare surface of sec- 
tion plot, it is in itself a difficult problem in more than 
two degrees of freedom. The geometric appearance of 
the orbits in configuration space does not provide a use- 
ful guide. A major part of this Letter will therefore be 
devoted to the development of a viable method to assign 
winding numbers to individual POs. 

Once winding numbers (and the associated action vari- 
ables, the calculation of which is then straightforward) 
are available, the situation in Fig. [5] can be character- 
ized by the following quantitative criteria: (i) The sta- 
bility angles cf)f of the FPOs S' ± (i.e. the phase angles 
of the unimodular eigenvalues of their stability matrices) 
describe the rotation that each FPO imposes upon its 
neighborhood. The winding ratios of the 2-torus POs 
converge toward (j)f /2n as the FPOs are approached, 
(ii) In the same limit, the action variable correspond- 
ing to the degree of freedom along the FPO converges 
to the action of the FPO, and (hi) the action variable 
for the motion transverse to the FPO, which is given 
by the area the (original) invariant torus encloses in the 
Poincare plane, tends to zero. We will use these three 
criteria to identify analogous situations in higher dimen- 
sions, and we will show that in a similar manner families 
of 2-torus POs can themselves serve as organizing cen- 
ters for families of 3-torus POs. Notice that although 
the choice of action-angle coordinates and the associated 
winding numbers is not unique, within the hierarchical 
structure the organizing center imposes a specific coor- 
dinate system (along/transverse to the FPO) upon the 
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family it organizes. 

POs that are remnants of the same torus have nearly 
identical periods and actions, thus each point in Fig. 
represents POs generated from a single torus. As in the 
two-dimensional case, 2-torus POs occur in doublets. By 
contrast, 3-torus POs arise in groups of four 0|. At the 
basis of the hierarchy of POs in the crossed-fields hydro- 
gen atom at E = — 1.5 a.u. there are the FPOs for 1-torus 
POs) S + and S~ that were described in 0,0 ■ Tnc y 
are shaped nearly like Keplerian ellipses in the invariant 
x-y-plane. 2-torus and 3-torus POs are arranged in sc- 
ries with positive or negative slopes in the T-S'-diagram 
of Fig. ^ which form the T£ 3 and T-f 3 families, respec- 
tively. The POs contained in the x-y plane and illus- 
trated in Fig. Inform the family T$. The series of 2-torus 
POs end at boundary lines marked by the FPOs and 
their repetitions. This observation suggests that each of 
the FPOs serves as an organizing center for both families 
T% and T£- We will show below that this is indeed the 
case. (Notice that in three degrees of freedom each FPO 
has two transverse degrees of freedom characterized by 
two different stability angles. It can therefore organize 
two different families of 2-torus POs.) In addition, we 
will see that the families T^' 11 themselves act as orga- 
nizing centers for the families T£' n of 3-torus POs. The 
entire hierarchy thus obtained is illustrated schematically 
in Fig. At the energy E = —1.4 a.u., above the ion- 
ization saddle point, we find the same general structure, 
except that the FPO S~ is surrounded by a region of ion- 
izing trajectories and thus cannot serve as an organizing 
center for the Xlf' 11 . 

Winding numbers for the 2-torus POs and TJf can 
be assigned as follows: Counting the series in Fig.^yields 
the first winding number w\ . The second winding num- 
ber W2 is found by counting the POs within one series 
from bottom to top. It remains only to determine with 
what values the counting is to start in both cases. To 
this end, we use Fourier expansions of the coordinate 
functions x(t) and z(t). These spectra show two major 
peaks. One of them agrees with the series number w\ 
(starting with i«i = 1 in the leftmost series), whereas 
the location of the second is identified with the second 
winding number W2 (starting with W2 = 2 for the short- 
est orbit in the leftmost series). The Fourier expansions 
thus confirm the validity of the simple counting scheme 
described above. 

For the 2-torus POs we obtain the action coordinates 
I shown in Fig. 0] together with the actions and stabil- 
ity angles 4>i.2 of the FPOs S ,± . The limiting values for 
high winding ratios of the and families coincide 
with 4>^[ /2tt and <^~/27r of S~ , respectively. As the max- 
imum winding ratio is approached, I2 converges towards 
the action of the FPO, whereas I\ vanishes. According 
to the three criteria listed above, we can therefore con- 
clude that the FPO S~ serves as an organizing center for 
the two families and T% of 2-torus POs, and we can 
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FIG. 4: (color online). Action variables for the T% (red 
plus symbols) and the T£ (green triangles) families at E = 
— 1.5 a.u., F = 0.5 a.u. Black dots show the action and the sta- 
bility angles of the organizing center S~ . The inset presents 
the same situation in a different action-angle coordinate sys- 
tem, in which S + can be identified as another organizing cen- 
ter of the T 2 p ' n families. 

identify the angles conjugate to I2 and I\ as the degrees 
of freedom along and transverse to the central FPO S~ . 

At the lower limits of the winding ratios in Fig. 0] none 
of the action variables converges to zero, and it seems, 
therefore, that, contrary to the hypothesis derived from 
Fig. ^ the T% and families are not organized by S + 
the way they are by S~ . However, this difference is an 
artifact of the action-angle coordinate system. In fact, 
Fourier spectra suggest that for 2-torus POs close to S + 
a set of winding numbers w' should be used that starts 
the counting scheme with W2 = 4 for the shortest orbit 
in the leftmost series. This set of winding numbers leads 
to the action variables I' used in the inset of Fig. In 
this coordinate system, the approach to S + satisfies the 
three criteria set up above, whereas the approach to S~ 
does not. Therefore, both FPOs S + and S~ serve as 
organizing centers for the two families of 2-tori POs in 
an equal manner. (The POs in Fig. ^are thus labeled in 
the coordinate system induced by the nearest FPO.) 

For the TJf , the conclusion that S + and S~ serve as 
organizing centers for families of 2-torus POs has already 
been reached from the Poincare surface of section plot in 
Fig. [3 Notice that in spite of its intuitive appeal the 
Poincare plot is somewhat misleading because it seems 
to show two distinct families of 2-torus POs. In fact, 
however, the T-S'-plot of Fig. ^ indicates, and the as- 
signment of winding numbers confirms, that all planar 2- 
torus POs form the single family Xf. A two-dimensional 
Poincare surface of section is entirely unsuitable to an- 
alyze the geometry of the T% because these POs do not 
lie in a two-dimensional plane. Nevertheless, the quan- 
titative method developed here demonstrates that the 
relation of Tjf to the FPOs S ,± is the same as that of the 
E\ '. Apart from avoiding this ambiguity, our method for 
the analysis of phase space topology is also more general 
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FIG. 5: (color online), (a) Winding ratios for POs from the 
T p families at E — — 1.4a.u., F — 0.5 a. u. with the same 
symbols as in Fig. The most prominent resonance lines 
[EmiH are indicated. The inset shows POs from the T n 
families, (b) Action variable Is for the family. For clarity 
the dotted line 7 from (a) has been subtracted on the vertical 
axis. 

than Poincare surface of section techniques in that it can 
now readily be applied to the 3-torus POs. 

The 3-torus POs in Fig. 2] fall into the same series 
as the 2-torus POs, but they possess a third winding 
number 103. The latter manifests itself in an additional 
major peak in the Fourier spectra of the coordinates, and 
therefore can be assigned by a straightforward extension 
of the technique used to classify the 2-torus POs. 

Fig. |S{a) shows the winding ratios of 2- and 3-torus 
POs. For the 2-torus POs, the missing ratio W3/W2 is re- 
placed with the stability angle that describes the dynam- 
ics transverse to the original torus. The stability angles 
of the 2-torus POs limit the winding ratios of the 3-torus 
POs. As this limit is approached, the action variable I3 
of the 3-torus POs tends to zero, as shown in Fig.OJb), 
whereas I± and I2 converge toward the values obtained 
from the 2-torus POs. These three observations are en- 
tirely analogous to the three criteria laid out above that 
characterize the relation of a family of 2-torus POs to its 
organizing FPOs. We can therefore conclude that we are 
here presented with a higher-dimensional generalization 
of that situation: The T2 families serve as organizing cen- 
ters for the 3-torus POs, and they impose a distinguished 
coordinate system onto the T 3 in which W3 and the action 
^3 identify the direction transverse to the Ti- 

In summary, we have demonstrated how periodic or- 



bits in a non-integrable multidimensional Hamiltonian 
system can be used to reconstruct a hierarchy of phase 
space structures that is organized by a stable PO at its 
center. The POs on each level of this hierarchy serve as 
organizing centers for the POs of the next-higher level. 
We have derived three general criteria to diagnose this 
relation between POs. They arc independent of the un- 
derlying dynamical system and can therefore be expected 
to be of wide applicability. At the same time, the calcu- 
lation of the associated action variables provides the key 
prerequisite for an Einstein-Brillouin-Keller torus quan- 
tization that the crossed-fields hydrogen atom has so far 
resisted. It thereby paves the way to an immediate, and 
important, semiclassical application of the purely classi- 
cal results obtained here. These results will be presented 
in a forthcoming publication. 
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